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Abstract 

The real time evolution of non-equilibrium expectation values with soft length 
scales ~ k^ 1 > (eT) _1 is solved in hot scalar electrodynamics, with a view to- 
wards understanding relaxational phenomena in the QGP and the electroweak 
plasma. We find that the gauge invariant non-equilibrium expectation values 
relax via power laws to asymptotic amplitudes that are determined by the 
quasiparticle poles. The long time relaxational dynamics and relevant time 
scales are determined by the behaviour of the retarded self-energy not at the 
small frequencies, but at the Landau damping thresholds. This explains the 
presence of power laws and not of exponential decay. In the process we red- 
erive the HTL effective action using non- equilibrium field theory. Furthermore 
we obtain the influence functional, the Langevin equation and the fluctuation- 
dissipation theorem for the soft modes, identifying the correlators that emerge 
in the classical limit. We show that a Markovian approximation fails to de- 
scribe the dynamics both at short and long times. We also introduce a novel 
kinetic approach that goes beyond the standard Boltzmann equation by in- 
corporating off-shell processes and find that the distribution function for soft 
quasiparticles relaxes with a power law through Landau damping. We find 
an unusual dressing dynamics of bare particles and anomalous (logarithmic) 
relaxation of hard quasiparticles. 
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I. INTRODUCTION 



There is currently a great deal of interest in understanding non-perturbative real time 
dynamics in gauge theories at high temperature, both within the realm of heavy ion collisions 
and the study of the quark gluon plasma |]- |§, as well as the possibility for anomalous 
baryon number violation in the electroweak theory [p|JT0[] . In both situations the dynamics 
of soft gauge fields with typical length scales > (gT)^ 1 is non-perturbative. 

Their treatment requires a resummation scheme where one can consistently integrate 
out the hard scales associated with momenta ~ T to obtain an effective theory for the soft 
scales. This is the program of resummation of hard thermal loops [IT]- ]T5|] . Physically, the 



hard scale represents the typical energy of a particle in the plasma while the soft scale is 
associated with collective excitations |16| . 

The recognition of the non-perturbative physics associated with soft degrees of freedom 
has led to an effort to describe the dynamics by implementing numerical simulations of 
classical gauge theories [[Lj]]- since soft degrees of freedom have very large occupation 
numbers and could in principle be treated classically ||20|| . Effective classical descriptions 
for the infrared bosonic modes have been obtained consistently in scalar field theory by 



integrating out the hard modes p5| . However, it was recognized that the dynamics of the 
soft modes in gauge theories is sensitive to the hard modes P3-EU[ and that the Rayleigh- 
Jeans divergences associated with the hard modes provide non-trivial contributions to the 
soft dynamics. 

For example, the one-loop correction to the gauge boson self-energy contains a leading 
contribution from the hard momenta of order g 2 T 2 which gives the HTL (hard thermal loop) 
contribution, and a subleading contribution from the soft scales. The T 2 dependence is a 
reflection of the UV quadratic divergence of the zero temperature theory which is cutoff at 
momentum scales > T by the Bose-Einstein factor. Restoring the appropriate h's we see that 
this contribution is of order ft(g 2 T 2 / ft 2 ) where the first ft is associated with the loop and the 
denominator follows from the usual manner in which ft enters with temperature. This term 
is therefore 0(T 2 /ft) and reveals the usual Rayleigh- Jeans divergence. For hard external 
momenta (K ext ~ T) this one-loop correction to the propagator is obviously subleading and 
bare perturbation theory is valid. But when the external momenta are soft (K ext ~ gT), 
clearly the one-loop correction is of the same order as the tree level term. This is in fact at 
the heart of the breakdown of the perturbative expansion. The problem is resolved by using 
HTL-resummed propagators and vertices for the soft external lines while hard scales may 
always be treated within the usual perturbation theory. This procedure is akin to obtaining 
a Wilsonian effective action for the soft modes by integrating out all the momenta above a 
certain soft scale which in this case is gT. For a detailed discussion of the relevant issues we 



refer the interested reader to the original works of Braaten and Pisarski [11,12 



All of the HTL contributions may be divided into two distinct categories: i) the con- 
tributions from tadpole diagrams and ii) those from diagrams with discontinuities. While 
the tadpole contributions are independent of the external momenta, the diagrams with dis- 
continuities lead to momentum dependent terms and it is these that lead to the non-local 
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effective HTL Lagrangian. The non-locality of the HTL effective Lagrangian for the soft 
modes originates in the process of Landau damping which results in discontinuities below 
the light cone [ITJ- [ p^ , p0| -|32[| . A very interesting method to deal with the non-locality in 



the HTL effective action in numerical simulations of classical gauge fields has been recently 
proposed [|33| and is based on the particle method akin to that used in transport theory. 



This method has been used to study the diffusion of Chern-Simons number in a lattice ap- 
proach ||34|| . A local Hamiltonian approach that is intrinsically gauge invariant has also been 



recently proposed an d has the potential for numerical implementation. A proposal to 
study the classical dynamics of soft gauge fields in terms of an effective Langevin equation 
has been put forth in P§| . However in our view such a proposal does not seem to incorpo- 



rate consistently the non-Markovian nature of the noise that is a result of the non-localities 
associated with Landau damping. 

The focus of this article is precisely to study in detail the real-time dynamics of the 
evolution of gauge fields, given an arbitrary field condensate in the initial state. This study 
will reveal that Landau damping processes dominate the most relevant aspects of the 
dynamics and as argued above, determine the non-local aspects of the HTL effective action. 
The main goal of this investigation is: i) to provide a deeper understanding of the time scales 
associated with dissipative off-shell processes, ii) a consistent microscopic description that 
can be used as a yardstick to test lattice results on real-time correlation functions, and iii) a 
detailed real-time description of relaxation and kinetics of soft collective excitations in gauge 
theories. 

Landau damping |[7|,p!T|-[T3| occurs when a hard quasiparticle from the thermal bath (with 



momentum ~ T) scatters off a soft collective mode (momentum ~ gT), borrowing energy 
from (and damping) the soft excitations in the process. Simple kinematics dictates that these 
processes can occur only off shell and below the light-cone i.e. for spacelike four-momentum. 
Furthermore, Landau damping gives a non-zero contribution only in the presence of a heat 
bath and when the external momentum is non-zero. Phrased differently, the Landau discon- 
tinuities are purely thermal cuts arising only at non-zero temperature and lead to damping 
of spatially inhomogeneous field configurations only. Whereas the real-time dynamics of sim- 
ilar processes has been studied in a scalar theory [Q, such a study is lacking for the case 



of gauge fields. The gauge boson self-energy in the presence of these processes has been 
known for a long time ppj-p 7 ^ and can be computed in the imaginary time formalism of 
finite temperature field theory in the HTL limit pTfl . 

The main focus and goals of these article are: 

• To compute explicitly the real time evolution of inhomogeneous field configurations 
(non-equilibrium field expectation values) in the ultrarelativistic plasma as an ini- 
tial value problem. We linearize the field equations of motion in the condensate 
amplitude. In this weak field regime, the evolution equations for the condensate can 
be solved in closed form through Laplace transform. The analytic structure of the 
propagator in momentum space (s-plane) determines the real time behaviour of the 
solution. We obtain the propagator to one-loop order in the HTL approximation. In 
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this approximation the long-time behaviour of the condensate turns out to be governed 
by the Landau discontinuities resulting in Landau damping processes. Although the 
HTL corrections to the gauge boson self-energy are well-known, we believe that the 
calculation of the real time dependence of the damping of soft excitations is new. 

While an understanding of the real-time relaxation of non-equilibrium, inhomogeneous 
field configurations is of fundamental importance in the physics of relaxation in the 
QGP 0, such a calculation also has phenomenological implications for sphaleron in- 
duced B-violating processes. It has been recently pointed out [f27|-|29| that standard 
estimates of the topological transition rate at finite temperature in the electroweak the- 
ory ignore the effects of damping in the thermal bath and these authors have argued 
that Landau damping plays a very important role. Since the sphaleron is an inhomoge- 
neous excitation associated with a soft length scale (~ l/g 2 T) Landau damping effects 
and the full HTL-resummed propagators must necessarily be taken into account when 
studying the sphaleron damping rates p7|-p^ . 



One of our main goals is to assess in detail the real time non-equilibrium dynamics 
of soft excitations in the plasma with particular attention to a critical analysis of the 
long-standing belief that the small frequency region of the spectral function dominates 
the long time relaxational dynamics. We find, to the contrary that the Landau damp- 
ing thresholds at uj = ±k determine the long- time dynamics and that the early time 
dynamics is sensitive to several moments of the total spectral density. This is an im- 
portant point that bears on recent arguments that seek to clarify the damping effects 



on the sphaleron rate [|27H29|j. We analyze this novel result in detail both analytically 
and numerically, thus proving that the long time behavior is dominated by the Lan- 
dau damping thresholds and that the small frequency region gives rise to sub-leading 
corrections to the long-time dynamics in the leading order HTL approximation. 

In this article we concentrate on the case of scalar electrodynamics (SQED) since this 
theory has the same HTL structure (to lowest order) as the non-abelian case p7|p8 



Most of our results can therefore, be taken over to the non-Abelian case with little or 
no changes at least in the lowest order HTL approximation. Scalar electrodynamics 
has already been used as an example to study the HTL resummation for the infrared 
modes f38| , but the scope of this article is different in that we study explicitly the real 



time dynamics of the damping processes. 

After resummation of the one-loop HTL contributions, at long times the relaxation of 
either transverse or longitudinal field expectation values is given by two contributions. 
The first is from the quasiparticle modes, and is standard - an oscillatory function in 
time. The second contribution arises from branch point singularities in the HTL self 
energies at nonzero frequency. These produce correlations in time which are oscillatory 
times power law tails; these power law tails are a new feature of HTL's. 

Long time power law tails in current-current correlators have been recently reported 



in reference 28 
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Furthermore we obtain consistently the effective Langevin equation for soft modes by 
integrating out the scalar fields and obtaining the influence functional for the gauge 
invariant observables. This allows us to extract the noise correlation function that dis- 
plays all of the non-localities associated with HTLs and Landau damping. By deriving 
the relevant fluctuation-dissipation relation we identify the proper correlation function 
that emerges in the classical limit. This analysis reveals that the range of both the 
dissipation kernel as well as the noise correlation function are determined by the soft 
scale. This results in the kernels being long-ranged, typically falling off with a power 
of time and with no Markovian limit. Again this result is deeply related to the non- 
localities of the HTL effective action and prevents a local description of relaxational 
dynamics associated with Landau damping. 

Our detailed analysis and the consistent formulation of the influence functional starting 
from the microscopic Lagrangian, unequivocally leads us to the conclusion that an 
effective stochastic Langevin description of gauge field relaxation to leading order in 
the HTL limit is non- Markovian. Because of the long-range kernels associated with 
Landau damping, a Markovian limit cannot be consistently extracted. This result 
points out the limited utility of a Langevin equation for describing relaxation via 
Landau damping. 

Having established the relaxation of inhomogenous gauge field configurations via off- 
shell processes associated with Landau damping, we ask how these processes contribute 
to the relaxation of the distribution function of transverse degrees of freedom. Clearly 
the evolution of the distribution function cannot be described to this order by a Boltz- 
mann equation, since this kinetic approach only includes on-shell processes. Thus we 
provide one of the novel results of this work: we incorporate the non-equilibrium relax- 
ation effects from Landau damping into a kinetic equation that describes the relaxation 
of the occupation number of transverse gauge fields. This kinetic equation incorporates 
off-shell effects and therefore constitutes an advance over the usual Boltzmann kinetic 
description in terms of completed collisions. We argue that since Landau damping 
results in an exchange of energy (and momentum) between the quasiparticles in the 
bath and the out-of-equilibrium field configuration, this in turn will naturally lead to 
a depletion of the particle number from the field configuration and must necessarily 
be included in any accurate kinetic description which aims to probe relaxational phe- 
nomena on the relevant time scales. A framework to study these transient, off-shell 



relaxational phenomena as initial value problems is given in refs. [j40|-|42l1. The anal- 
ysis presented in this article implies a Dyson-like resummation that goes far beyond 
Boltzmann kinetics. 

We compare the results from these new kinetic equations to several different approx- 
imations to the kinetics. In comparing the relaxation of dressed soft quasiparticles 
and that of bare particles and hard quasiparticles, we find a remarkable dressing dy- 
namics of the degrees of freedom in the medium, and anomalous relaxation for the 
hard quasiparticles. We also argue that a proper kinetic description must incorporate 
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consistently the HTL effects and the quasiparticle nature of the excitations. 

In Section II we introduce the model under study and focus on a gauge invariant de- 
scription of the dynamics. Section III is devoted to a derivation of the non-equilibrium 
equation of motion for the inhomogeneous condensate to one-loop order in the hard thermal 
approximation. The real time relaxational dynamics of the inhomogeneous configuration 
via Landau damping is investigated in detail in Section IV. Contact is established with the 
fluctuation-dissipation theorem and stochastic dynamics in Section V, where we derive a 
Langevin description for the soft gauge invariant degrees of freedom in the thermal bath 
and recognize the relevant correlation functions that emerge in a (semi) classical stochastic 
description. We introduce a new kinetic description of transport phenomena induced by the 
non-collisional Landau damping process in Section VI. In this section we study the relaxation 
of the distribution function for soft quasiparticles, bare particles and mention some interest- 
ing features of the relaxation for hard quasiparticles. Finally, we summarize our analysis, 
discuss the modifications that will arise when higher order corrections leading to collisional 
lifetimes are included and present our conclusions and possible future directions of study. 

II. PRELIMINARIES 

As mentioned in the Introduction, our ultimate goal is to understand relaxational pro- 
cesses associated with off-shell effects, such as Landau damping, in a Non-Abelian gauge 
theory [[7J]. What we will do here is to treat the same problem in the context of scalar 
quantum electrodynamics (SQED) model. To leading order, we expect that this should be 
a good analogue of what happens in the Non-Abelian case; it also has the advantage that 
it is simpler to deal with, and as we will see below, it can be cast from the outset in terms 
of gauge invariant variables. This will eliminate any ambiguities associated with the usual 
problem of gauge dependence of off-shell quantities. 

We will start with some inhomogeneous field configuration which is excited in the SQED 
plasma at t — 0. What we want to do is to follow the time development of this configuration 
as it interacts with the hard modes in the plasma, and in particular, we want to know 
whether the relaxation is of the usually assumed exponential sort or something different. 

Let us first reformulate SQED in terms of gauge invariant variables. The SQED La- 
grangian is given by 

£ = - m 2 |$| 2 - -F^F^ , 

= (0„ - ieAJ® . (2.1) 

A description of the dynamics in terms of gauge invariant observables begins with the 
identification of the constraints associated with gauge invariance. The Abelian gauge theory 
has two first class constraints, namely Gauss's law and vanishing canonical momentum for 
Aq. What we will do is project the theory directly onto the physical Hilbert space, defined as 
usual as the set of states annihilated by the constraints. The procedure is simple. First, we 
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obtain gauge invariant observables that commute with the first class constraints and write 
the Hamiltonian in terms of these. All of the matrix elements between gauge invariant states 
(annihilated by first class constraints) are the same as those that would be obtained by fixing 
Coulomb gauge V • A = 0. The Hamiltonian, when defined in the physical subspace, can be 
written solely in terms of transverse components and includes the instantaneous Coulomb 
interaction as would be obtained in Coulomb gauge. This instantaneous Coulomb interaction 
can then be traded for a gauge invariantLa.gra.iage multiplier field A (x, t) (a non-propagating 
field whose canonical momentum is absent from the Hamiltonian) linearly coupled to the 
charge density p(x,t) and obeying the algebraic equation of motion V 2 A (x, t) = p(x,t). 
Alternatively, one can use a phase-space path integral representation of the generating func- 
tionals, trade the Coulomb interaction with a Lagrange multiplier linearly coupled to the 
charge density and perform the path integral over the canonical momenta as usual. Both 
methods lead to the following Lagrangian density: 

e 2 A T ■ A T $ f $ + 

(2.2) 

where At is the transverse component of the gauge field. 

In order to provide an initial value problem for studying the relaxational dynamics of 
charge density fluctuations we introduce an external source J7l(x, £) linearly coupled to A 
and study the linear response to this perturbation. Furthermore, it is convenient to intro- 
duce external sources coupled to the transverse gauge fields to study the linear response of 
transverse gauge field configurations. These external fields could in principle play the role of 
a semiclassical configuration coupled to small perturbations in a linearized approximation. 
Therefore we include external source terms in the Lagrangian density, 

C -> C - J L {x, t)A (x, t) - J T {x, t) ■ A T (x, t) . (2.3) 

The relaxational dynamics of our initial inhomogeneous configurations is clearly an out 
of equilibrium process, and needs to be treated by an appropriate formalism p3|-|45[] . 

In the Schrodinger picture the dynamics is completely described by a time-evolved density 
matrix p that obeys the quantum Liouville equation: 

ih^ = \p,H], (2.4) 

where H is the Hamiltonian of the system. The expectation value of any operator O is given 
by 

(O) = Tr[p(t)0] . (2.5) 

For thermal initial conditions with an initial temperature given by 1/(3, the density matrix 
at t = is pi — e~P Hi and the above expectation value can be rewritten easily as a functional 



£ = + ~d„A T ■ d»A T - eA T ■ j T 

+i (VA,) 2 + e 2 A 2 Q $t$ _ Z eA) 
j T = i($ty r $ - V T $ f $) • 
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integral defined on a complex-time contour. Notice that Hi is not the Hamiltonian of the 
system for t > 0. The system thus evolves out of equilibrium. The contour has two branches 
running forward and backward in time and a third leg along the imaginary axis stretching 
to t = —if3. This is the standard Schwinger-Keldysh closed time path formulation of non- 
equilibrium field theory (see ref. pEBj-pEBI] for details). Fields defined on the forward and 
backward time contours are accompanied with (+) and (— ) superscripts respectively and 
are to be treated independently. The expectation value of any string of field operators may 
be obtained by introducing independent sources on the forward and backward time contours 
and taking functional derivatives of the generating functional with respect to these sources. 
The imaginary time leg of the complex time contour does not contribute to the dynamics. 
Since the path integral represents a trace, the initial and final states must be identified 
and therefore all the local bosonic fields 0(x, t) satisfy the Kubo-Martin-Schwinger (KMS) 
periodicity condition 



The non-equilibrium SQED Lagrangian is given by 



C 



noneq 



c 



c 



(2.6) 



(2.7) 



Perturbative calculations are carried out with the following non-equilibrium Green's func- 
tions: 

• Scalar Propagators: 



where (a, b) €{+,—}. 

g++(m') = G>(t,t')e(t-t') +G<(t, t')e(t'-t) , 

GjTftO = G>(t,t')e(t'-t)+G<(t,t')Q(t-t') , 
Gf*(f,0 = -G<>\t,t') , 

G>{t, t') = -i- [(1 + n k ) e -^(*-*') + n k e ^*(*-0 

G<{t,t') = ^- k e-^-V + (l + n fe ) e*^*-*) 
2uj k L 

1 



(2.8) 

(2.9) 
(2.10) 

(2.11) 



LU k = V k 2 + m 2 ; n k 



Photon Propagators: 



{AS(x,t)4}(x,t')) 



—i 



r1 3 Jc - 
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g++(k;t,t') = p ij (k) [g>(t,t')G(t-t f ) + g<(t,t f )e(t f -t)\ , 

grr^tX) = v VJ {k) [g>(t,t')e(t' -t) + g<(t,t')e(t-t')} , 
g^(k;t,t') = -v lJ (k) £< (>) (m'), 

g£(t,t') = ^ I (1 + N ^ e ~ tHt ' t ' ) + N k e lk{t -^] , (2.12) 
#(t,f ) = ^ e ~ i/£( *" 4,) + ^ + ^ ' ( 2 - 13 ) 



Here Vij(k) is the transverse projection operator: 

TV*) = ^ - ^ ■ (2.14) 
With these tools we are ready to begin our analysis of non-equilibrium SQED. 

III. LINEAR RELAXATION 

We 

now introduce the inhomogeneous non-equilibrium expectation values AT{x,t) ; Ao(x,t) 
which are excited at time t = in the plasma. The non-equilibrium expectation values of 
the transverse components represent electric and magnetic fields, whereas the expectation 
of the Lagrange multiplier field A$ corresponds to an initial charge density in the system. 

The dynamics of these non-equilibrium expectation value will be analyzed by treating 
At(x, t); Ao(x, t) as background fields, i.e. the expectation values of the corresponding fields 
in the non-equilibrium density matrix, and expanding the Lagrangian about this configura- 
tion. Therefore we split the full quantum fields into c-number expectation values (which are 
the non-equilibrium expectation value) and quantum fluctuations about these expectation 
values: 

A^\x , t) = A T (x, t) + 5A { t\x, t) ■ A^\x, t) = A (x, t) + SA^ } (x, t) (3.1) 
A T (x,t) = (A^\x,t)) ; A (x,t) = (A^\x,t)) (3.2) 

where the expectation values of the field operators are taken in the time-evolved density 
matrix. 

The equations of motion for the background field can be obtained to any order in the 
perturbative expansion by imposing the requirement that the expectation value of the quan- 
tum fluctuations in the time evolved density matrix vanishes identically. This is referred to 
as the tadpole equation [|5[] which follows from Eq.( [3.1|) and Eq.( [3.2|) : 



(Sffl) = ; (SA^) = . (3.3) 
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The equations obtained via this procedure are the equations of motion obtained by variations 
of the non-equilibrium effective action. The perturbative expansion needed to compute the 
relevant expectation values is obtained by treating all the linear terms [ |5| in the fluctuations 
as interactions along with the usual interaction vertices. 

Although in principle the tadpole method could be used to study arbitrary background 
configurations including non-perturbative ones (for e.g. sphalerons in the non-Abelian case) 
we restrict our discussions to the small amplitude regime in close analogy to the work of ref. 
|36| for scalar field theory. In other words, the effective action equation of motion will be 
studied in the linear approximation for the condensate amplitude so that O (A^) and O (Al) 
and higher orders will be neglected. In addition, the evolution kernel will be approximated 
to the one-loop order only. 

Defining the Fourier components of the electromagnetic condensate as 



A Ti (k, t)= I d 3 x e lk ' s A Ti {x, t), (3.4) 
Ao(k, t)= I d 3 x e its A (x, t) (3.5) 



we obtain the following equation of motion to one-loop order for the transverse part: 
d 2 

—A T i(k } t) + k 2 A T i(k,t) + 2e 2 {&<S>} A Ti (k,t) + 

-2e 2 dr — — p Ti p Tj [(1 + n p + n p+k ) sm{(u k+p + u p )(t - r)} 

JO J {Z7T) U} p LO p+ k 

+ {n p - n k+p ) sin{(u; fc+ p - u p )(t - r)}] A T j{k,r) = J T i{k,t) . (3.6) 

Here pt% refers to the component of the spatial momentum p which is transverse to the 
wave- vector k and Arj{k, r) stands for the Fourier transform of the external source ATj{x, t) 
in Eq.Q. 

The tadpole term which appears due to the 4-point 'seagull' vertex can be evaluated 
easily in the limit where T » m and yields the following hard thermal loop contribution 

, . , . n f 1 + 2n k d 3 k e 2 T 2 . . 

We remark that in evaluating the tadpole diagram above and in all subsequent calcula- 
tions it is always assumed that the zero temperature divergences have already been absorbed 
in the proper mass and wavefunction renormalizations. 

The longitudinal part obeys a similar equation: 



k 2 A (k,t) + e 2 J* dr J 



d 3 p 

o J {2^f 



"""" P - 1 ] (1 + n p + n p+k ) sin{(u; fe+ p + u p ){t - r)} 



+ ( + l)(n p - n k+p ) sin{(w fc+p - uj p )(t - r)} 



UJ p 



A (k,r) = J L (k,t) . (3i 
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It should be noted that the equation of motion for the longitudinal component has no 
time derivatives, indicating the non-dynamical nature of the field. The source term for 
the longitudinal component is interpreted as an external disturbance that induces a charge 
density fluctuation in the SQED plasma. The response to a general disturbance can be 
obtained by convolution from the result of linear response to the impulsive perturbation. 

We also note that unlike the transverse case, there is no contribution from the tadpole 
term to the effective action equations of motion. This fact will be important in understanding 
the origin of a non-zero Debye mass. 

The nonlocal terms in Eqs. (|3.6|) and ( |3.8|) are the one-loop self-energies for the transverse 
and longitudinal components respectively, resulting from the photon-Higgs trilinear coupling. 
The first piece in the nonlocal terms proportional to 1 + n p + n k+p is the difference of the 
following creation and annihilation processes in the medium |J5| : 7 — ► $ $ with Bose 
factor (1 + n p ) (1 + n p+ k) and <5 $ — > 7 with a statistical factor 
proportional to 



Tl r) Tlr 



The piece 



(n p — rik+p) is the Landau damping contribution [jn^fKJ. When k = this 
contribution vanishes indicating that it affects inhomogeneous excitations only. Furthermore, 
it has no zero temperature counterpart since the Bose factor n k — n k+p vanishes identically 
at T = 0. This term arises from the difference of the processes 7 $ — > $ with statistical 
factor (1 + rik+p) n p and $ — > 7 $ with the factor (1 + n p ) n k + p . Although it does not give 
rise to an imaginary part for the on-shell self-energy, it will nevertheless have an effect on 
the physical processes associated with the relaxation of the inhomogeneous condensate as 
described in detail below. 

We can solve Eqs.( 
formed fields 



d[ |3.8|) via the Laplace transform. Introducing the Laplace trans- 



A T i(k,s) 
Ao{k, s) 



dt e- st A Ti (k } t) 



dt e- st A (k,t) 



(3.9) 
(3.10) 



and performing the transform on the above equations of motion we get following the same 
methods as in ref. |3l) , 



{s 2 + k 2 + e 2 T 2 /6)A Ti (k,s) + 
2 / d 3 p 



-2e z 



+ (n p - n k+p ) 



(2ir) 3 uJpUJ k+ p 



1 + n p + n p+k ) 



^k+p + Up 



{^k+p 



LO r , 



(3.11) 



s 2 + (uj k+p + Up) 2 



PTi Ptj A~Tj(k, s) = sA T i(k, 0) + A T i(k, 0) + Jn{k, s) 



for the transverse part, and 



k 2 + e 2 



d 3 p 



^k+p + Up 



0J„ 



s 2 + {u k+p + u p 



(3.12) 
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h 1 n p - n fc+p J 

cj p j s 2 + [uj k+p - u p y 



A (k,s) = J L {k,s) 



for the longitudinal component. The Laplace transform variable s plays the role of an 
(imaginary) time component of the photon four-momentum. 

IV. REAL TIME LANDAU DAMPING 

With the equations of motion for the non-equilibrium expectation value in hand, we turn 
to the analysis of its time evolution, paying particular attention to those parts of the spectral 
density which contribute to the long time dynamics. 

We do this by solving Eqs. (|3.11| , |3.12|) via the inverse Laplace transform. Recall that 



this requires an integration along a contour parallel to the imaginary s axis placed in such a 
way so as to have all the singularities of the integrand to the left of the contour. To do this, 
we will need a clear understanding of the analytic structure of the photon self-energy in the 
s-plane. We will explicitly outline the details of the calculation for the transverse part only, 
since those for the longitudinal part are similar. 

A. The Transverse Part: 



It is illuminating to write the equation of motion (|3.11| ) in terms of a spectral density 
function (u/, k) [p^] , 



-> r+oo tjj' _ 

(s 2 + k 2 + e 2 T 2 /6)A Ti (k,s) + du' Pij{u',k) A Tj (k,s) (4.1) 

J— oo S -T UJ 

sA Ti (k, 0) + A Ti (k, 0) + J Ti (k, s) 



where 



Pij (u', k) = -2e 2 [ d V p Tl p Tj x (4.2) 

x [(1 + n p + n p+k )5(uj' - uj k+p - uj p ) + (n p - n k+p )5(uj k+p -u p — to')} . 
The one-loop transverse self-energy in the s-plane is given by: 

e 2 T 2 r°° lo' 

Vijik^is) = — Vij(k) + / dw'^— — Pij (u',k,s), (4.3) 

O J —oo S -j- UJ 

where we have recognized the fact that pij oc Vij(k). 

Using Vij(k)ATj{k, s) = Ari{k,s) the Laplace transform of the transverse part of the 
condensate is given by: 
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The real time dependence of the inhomogeneous background is given by the inverse Laplace 
transform of the above expression which is in fact the retarded propagator defined in the 
s-plane. The inverse transform is calculated by performing the following integral along the 
Bromwich contour. 

c + io ° ds 



_, rc+ioo _, 

A Ti (k,t) = —e st A Ti (k,s 

J c—ioo Zlll 



c +*°° ds_ st s A Ti (k, 0) + A Ti (k, 0) + J Ti (k, s) 
-too 2m 6 s 2 + k 2 + E*(s) ' 1 ' ' 

Here we choose c > 0, such that the contour is to the right of all the singularities in the 
s-plane. 

The time dependence of this integral crucially depends on the analytic properties of the 
propagator and hence a clear understanding of the poles and cuts of the retarded propagator 
is essential. At this stage we will set the external current to zero, and analyze the contribution 
of the sources at the end of this subsection. 

In all of the above and in what follows it is implicitly assumed that the zero temperature 
divergences have been dealt with already by renormalizing the amplitude of the field, i.e. 
wave-function renormalization and that we are working with the subtracted spectral function. 
Recall that the divergences are determined solely by zero temperature fluctuations. 

We define the real and imaginary parts of the self-energy near the imaginary axis through 

± 0+) = E^O'u; ± + ) + iS* (iu ± + ). (4.6) 
In fact it is easy to see from Eq. Q4.3p that 



S*(io; ± 0+) = Tsgn(u)^( Pij (\Lu\) - Pij (-\u\)) (4.7) 



G f d p p Sill u 
±— — sgn(a;) / [(1 + n p + n k+p )S(u; k+p + u p - H) + 

07T J UJ p UJk+p 

+ (n p - n k+p ) {5{u k+p -u p - \u\) - 5{u k+p -u p + . 



This shows that ^(iu) — + ) = — ^(iu + + ) and that the cuts will appear whenever the 
delta functions are satisfied. Analysis of the arguments of the delta functions reveals that 
the two-particle cuts arising from the first term in Eq.(I77) stretch from s = ±i(m + uj k ) to 



s = ±ioo. The second and third terms in Eq. (|4.7|) contain the hard thermal contributions. 
They have support when u> k+p — u p ~ k ■ p = ±u) which corresponds to the four-vector 
(u, k) being spacelike. The resulting branch cut runs from s = —ik to s — +ik. Hence 
these processes are induced only by off-shell space-like photons. These HTL contributions 
lead to Landau damping. In summary, the physical region singularities and the Landau 
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discontinuities show up as discontinuities in the imaginary part of the self-energy S*(s) when 
approaching the imaginary axis of the s-plane i.e. s = iu. 

In addition, as argued in the previous section the contribution of the physical region cuts 
from s = ±i(m + u)k) to s = ±ioo to the self-energy are oc InT only, and therefore subleading 
compared with the 0(T 2 ) terms. 

Furthermore the long time dynamics due to the cuts will be dominated by the thresholds 
(or the end-points). Since the end-points of the Landau damping discontinuities are at s ~ 
±ik, they will be the dominant contribution and the two particle cuts from s = ±i(m+Uk) to 
s = ±00 will be subleading at long times. Thus we can simply focus on Landau damping both 
as the leading high temperature and long time contributions. This argument is necessary 
because it is not a-priori obvious that the high temperature and long time limits are described 
by the same processes. 

The leading term in the self-energy oc T 2 is given by, 



S*(io;±0 + ) = 
e 2 

87T 



(4- 



dp p 2 



dn p 
dp 



) I dx (1 — x 2 )kx [5(kx — \oj\) — 5(kx + \u\ 



-1 



The real part can be obtained either by using dispersion relations or by explicitly solving for 
the hard thermal self-energy by calculating the relevant integrals and we find 



e z T 
12 



s 2 ,s /, s 2 \, (is-k s 

2 ¥ + i k[ 1 + ¥ ln [iITk 



+ e>(inT). 



(4.9) 



From this expression the transverse self-energy along the imaginary axis when approaching 
from the right, can be read off easily: 




.e 2 T 2 7r to / u. t , : 



00 



(4.10) 



and agrees with known results [31 , 32 , 37 , 38 , 46 



In addition to the branch cut singularities, the retarded propagator will also have isolated 
poles corresponding to the quasiparticle excitations which can propagate in the plasma. The 
poles for the transverse excitations will be given by the solutions to 



fr 2 
12 

e 2 T 2 7r up 



uJ 



r — k + 
+i 



,t4 
' k 2 




k 2 



j In 


k + up 


k — up 





+ 



u. 



1-^ e { k 2 -u 2 P ) 



12 k 

They will be in the physical sheet provided the imaginary part vanishes at the pole. 



(4.11) 
(4.12) 
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Though the equation cannot be solved analytically, in the case of interest when the 
external momenta are extremely soft (for e.g. k ~ e 2 T <C eT) representing a small amplitude 
long wavelength field configuration, the approximate location of the poles is found to be 

s = ±iuj P ^±if 

The two-particle cuts were shown to run from s = ±i(m + ujk) to ±200 where m is the 
mass of the scalar. A consistent HTL resummation should also include the shift in the scalar 
masses, thus ensuring that to this order the quasiparticle pole is in the physical sheet. Higher 
order contributions will provide a collisional broadening to the pole. It is a noteworthy point 
that the quasiparticle poles are located beyond the Landau discontinuities which stretch 
from — ik to +ik in the s-plane, and below the two-particle threshold. 

In summary, to this order in the HTL approximation, the analytic structure of the re- 
tarded propagator in the high temperature limit features a Landau discontinuity running 
from — ik to +ik and quasiparticle poles at s = ±iu p (k). The two-particle cut contributions 
have been shown to give a subleading contribution both in temperature and in the long time 
dynamics. 

Using Eq.( ^.5| ) we can now invert the transform by deforming the contour and wrapping 
it around the poles and the cuts to pick up the corresponding residues and discontinuities 
repectively so that, 



A Tl (k, t) = A$ e (k, t) + A^(k, t), 



(4.13) 



The contributions from the quasiparticle poles add up to give a purely oscillatory behaviour 
in time. The residues at the poles give rise to a wave function renormalization with both 
T-dependent and T-independent contributions. The T-independent contribution contains 
the typical logarithmic divergence and has been absorbed in the usual zero temperature 
wave-function renormalization. We thus obtain, 



A^f(k,t) = Z*[T] 



An{k, 0) cos(cjpt) + A T i(k, 0) 



sin(a;pt) 



Up 



Z l [T] 



dH\iuj) 
duj 2 



-i 



co=u:p—eT/3 



(4.14) 
(4.15) 



Here Z t (T) is the temperature dependent wave function renormalization defined on-shell 
at the quasiparticle pole whose leading HTL contribution is obtained from the self-energy 
( |4.10| ). The continuum contribution is given by 



2 rk S^(iu; + + ) A T i(k, 0) uj cos(cut) + A T i{k, 0) sin(cut) 
7T Jo dU [lu 2 -k 2 - X^Mf + m(iu + 0+)f 



(4.16) 



Evaluating Eqs. (|4.13| ),( |4.14D and (|4.16|) at t = we obtain an important sum rule, 



2 r l 
Z\T\ + - / 

7T JO 



ujYt^iuj + 0+) 



[u 2 - k 2 - S^Mf + [£*(zu; + 0+)f 



(4.17) 
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This sum rule is a consequence of the canonical commutation relations but its content in 
the HTL approximation is that in the high temperature limit the wave function renormal- 
isation which is evaluated on-shell is completely determined by the Landau discontinuities 
which originate from strongly off-shell processes. A similar sum rule was also obtained in 



II | using different methods. 

The integral over the cut QL16Q cannot be evaluated in closed form but its long time 
asymptotics is dominated by the end-point contributions as can be understood from the 
following argument. The integral along the real uj axis from uj = to uj = k can be 
obtained by deforming the integral into the upper complex uj plane so that it runs along 
uj — iz ; < z < oo, then around an arc at infinity and back to the real axis along the line 
uj — k + iz ; < z < oo for the term oc e tuJt and similarly into the lower complex plane for 
the term cx e~ tuJt . A detailed analysis of the asymptotic behavior of the integral reveals that 
only the uj = k end-point contributes because the contributions from uj = vanish for large 
t faster than any negative power of t. This is a consequence of the regular behavior of the 
spectral density in the vicinity of uj = 0. 

We find the contribution from the uj = k end-point in the long time limit t » 1/ A; to be 
given by, 



1 + 



(4.18) 



^From the sum rule ( [4.16 ) and ([4.17 ) we find that the early time behavior is approximately 
given by 



A Ti {k, t) = A Tl (k, 0) [l - e 2 T¥A e + 0(t 4 )] + A Tl (k, 0) A ^f t [1 + 0(t 2 )} , (4.19) 

where A 6i0 are constants depending on the wave-function renormalization and moments of 
the spectral density. 

At long times, the dominant contributions are from the nearest singularities in the com- 
plex uj plane. This includes the usual contributions from the quasiparticle modes at uj = ±ujp, 
given by Eq. ( 4.14|) . In addition, there is also the contributions from the branch points at 



uj = ±k, Eq.( [4.18| ). There is no contribution from uj = 0, because the HTL self energy is 
regular about zero frequency. 

Figure (|1|) shows the cut contribution A^ ,t (k, t)/AT(k, 0) vs. time for e 2 T 2 /k 2 = 2 and 
figure (fj) shows t 2 x A^ (k , t) / Ar{k , 0) vs. t (in units of 1/k) for m 2 D /k 2 = 12 for the case 
A Ti (k,0) = 0. 

The point being made is that the real time dynamics of the condensate is completely 
determined by the analytic structure of the retarded propagator and the global structure of 
the spectral density in the s-plane. 

The second important point to note is that the long time behaviour is a power law ~ t~ 2 
(times oscillations) and not an exponential decay. This means that Landau damping 
effects cannot be reproduced by phenomenological 'viscous' terms of the type ~ T ^ neither 
at long nor at short times. The failure of such a phenomenologically motivated ansatz was 
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already noticed at zero temperature in different contexts in ref. ||44j| . We stress that such a 
description not only fails to reproduce the power law behaviour but in fact ignores all the 
non-local physics of Landau damping which is so clearly encoded in the Hard Thermal Loop 
kernels. 

One might argue that higher order processes, both Landau damping and collisional, could 
lead to an exponential relaxation. However, the point of the above analysis is to argue that 
the full relaxational physics will be described by a competition between the power laws from 
the lowest order Landau damping contributions and the higher order exponential damping. 
The time scale of interest will determine which process dominates. 

For non-zero external sources we can obtain the general real time evolution by inserting 
the Laplace transform of the source in Eq. (|4.5|) 



J Ti (k,s) = / dte~ st J Ti (k,t) . 
Jo 



This is an analytic function of s for Re(s) > provided Jri(k,t) is a non-singular function 
of time. Let us assume that Jn(k,s) is analytic for Re(s) > —a, where a > is a posi- 
tive number. [Since Jn{k,s) vanishes for Re(s) — > +00, Jri(k,s) must have singularities 
somewhere in the left half-plane. Otherwise, it will be an entire function which is zero at 
infinity and therefore identically zero]. The singularities of Sn(k, s) in the left-half s plane 
yield contributions to ATi{k,t) through Eq. ( |4.5| ) which decrease exponentially in time as 
e~ at . They are therefore subdominant compared with the Landau cut contributions (|4.18|) . 

The source term contribution from the small s = iuj region can be understood simply 
as follows. Within the above hypothesis, the source term Jn(k,iuj) can be expanded in a 
Taylor series for small uj 

J Ti (k, ioj) = J (k) + Jx{k) uj + 0(u 2 ) . 

The even term, Jo(k) can be absorbed into a shift of An{k,0), while the odd term 
Ji(k) can be absorbed into a shift of Ari{k,0), ( |4.5|) . Thus as before the asymptotic large 
time behavior is completely dominated by the branch points at uj = ±k. Under the above 
assumptions on the Laplace transform of the external current, we find the general time 
dependence to be given by 

Afffat)*^ (4.20) 

12 f ■ - 1 

- 2T2kf2 [[kAnik, 0) - Si(k)) cos(fct) + [A Ti {k, 0) + d{k)} sin(H)} 

Here Ci(u) and Si(u) stand for the Fourier cosine and sine transforms of the source J7r«(fc, t), 
respectively: 

roc roc 

Ci(uj) = dt cos(ujt) Jn(k,t) , Si(u>) = dt sm(ut) J Ti {k,t) . 
Jo Jo 
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We emphasize that to this order in the HTL resummation, the region near ~ of 
the spectral density is regular (no branch singularities) and therefore does not contribute to 
the long time dynamics, which is completely determined by the end point (branch point) at 
uj = k. 



B. Longitudinal part 

An analysis very similar to the one outlined above yields the following expressions for 
the longitudinal component of the photon self-energy in the case of an impulsive source 
i7l(x, t) = 5 3 (x)5(t) with spatial Fourier and Laplace transform given by JiXk, s) = 1. The 
case of a more complicated source can be obtained by convolution. In this case we obtain 



A Q (k,s) 



k 2 + E l {s) 



where 



e 2 T 2 



e 2 T 2 



6 ik 



In 



is 



is + k 



(4.21) 



(4.22) 



Thus the longitudinal self-energy along the imaginary axis in the s-plane, when approaching 
from the right is obtained as before to be: 



Z l (iuj + 0" 



e 2 T 2 



1 -In 

2k 



k + uj 



k 



UJ 



UJ 



(4.23) 



The location of the longitudinal quasiparticle or plasmon poles is given by the following 
dispersion relation, 



e 2 T 2 



k 2 + 



1 



UJ 

2k 



In 



k + uj 



k 



UJ 



.e 2 T 2 n uj 2 2 
i> — 7. — y<£>(k - uj ) = 0. 
6 k 



(4.24) 



For soft external momenta (k << eT) the plasmon poles can be seen to be at s = ±iuj ~ 
±ieT/3. The real time dependence of the longitudinal condensate is then found by inverting 
the transform using 



A (k,t) 



c + io ° ds 
—ioo 2ixi 



e st A (k : s) 



(4.25) 



where the contour is to the right of all the singularities as in Eq.( [4.5| ). As in the transverse 
case, at high temperatures the singularity structure is dominated by the discontinuity across 
the cut that runs from — ik to +ik corresponding to Landau damping, and the two plas- 
mon poles at ±iuj . The Bromwich contour is then deformed to pick up the cut and pole 
contributions so that 



A (k,t) = A p ole (k,t)+A^ ut (k,t), 



(4.26) 
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where 



A p ole (k,t) = -Z l [T] 



sm(u>ot) 



LO 



Z l [T] 



dY}{iu) 
duj 2 



(4.27) 
(4.28) 



u}=u) c^eT/3 



and 



A^(k,t) 



2 f k T}j{iuj + 0+) sin(utf) 

/ (lUl ty 

TT Jo ■ - 



k 2 + E l R (iuj) + E l j(icu + 0+) 



(4.29) 



Unlike the transverse components for which the sum rule is a consequence of the canoni- 
cal commutation relations, for the longitudinal component there is no equivalent sum rule 
because the field A (x) is a non-propagating Lagrange multiplier. 

The long time, asymptotic behaviour of the longitudinal condensate is dominated by 
the end-points of the integral Q4.29 ). The end-point uj = yields contributions that vanish 
for long t faster than any negative power of t, as it was the case for the transverse part 
( |4.16| ) . We find that the long time asymptotics for t > > 1/ k is dominated by the end-point 
contribution at uj = k, 



a 



cut 

asymp \ 



cut 

asymp 

12 



(k,t) 
a(x, t) = arctan 



(M) 



J^2 



dxe 



cos[kt + a(x, t)] 



1U & kt T 2 1U & kt T 16 



7T log 



*E1 _)_ Jog- 2 S 
, 



2 6Xp[2+ e^ 



(4.30) 

(4.31) 
(4.32) 



The integral in Eq. ( |4.30| ) cannot be expressed in terms of elementary functions and it 
is related to the v(x) function |4"9" |. For large t, one can derive an asymptotic expansion in 
inverse powers of log[/ct/c] by integrating Eq.( [4.30| ) by parts: 



cut 

''asymp 



(k,t) 



t— »oo 



12 cos(kt + (3(H)) 



e 2 TH 
(3(kt) = —arctan 



2 kt 



+ 



log 

TT log kt 

log 2 kt - 1 



1 + 



log* 



(4.33) 
(4.34) 



This expansion is not very good quantitatively unless t is very large. For example, 
for kt = 500, a>a% mp (k,t) given by Eq. (|4.30|) approximates AQ Ut (k,t) up to 0.1%, whereas, 



the dominant term in Eq. Q4.33D is about 30% smaller than a™t, mp (k, t). Figure (^|) shows 
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t[log 2 {kt/c) + 7c 2 /4}A^ ut {k,t) and (§ shows (12/e 2 T 2 ) cos[kt + /3(kt)], (see Eqs.(CT)-(^33|)), 
both figures should coincide for large enough t but we see numerically a sizable discrepancy 
even for t = 500 due to the slow convergence of the asymptotic expansion ( t4.33|) . 
A saddle point analysis yields an intermediate asymptotics of the form 



cut 

asymp 



(k,t) 



3\/7T cos(kt) 
e 2 T 2 t\n 2 - 5 [kt] 



(4.35) 



which although is a relatively good estimate within a time window, (see fig. (H)) is seen 
numerically to have a discrepancy of the same order as the dominant term above in a wider 
range of time. Nevertheless, the exact expression (|4.30|) is easily obtained numerically for 
arbitrary range of parameters. 

The case of a general source requires a convolution. The cut contribution takes now the 
form, 



7T JO 



k + 0+) [C l {uj) sin(cut) - S L {u) cos(ut)] 

dui 



k 2 + T} R {iuj) + Hftiu + O 



where Cl{oj) and Sl(uj) stand for the cosine and sine Fourier coefficients of the source 
J~L(k,t), respectively. 

For a longitudinal source fullfiling the same general analyticity assumptions as the trans- 
verse source we can expand in series for small uj as 

J L (k, iu) = JqL + J XL w + O(co 2 ) . 



As for Eq. ( |4.29| ) the end-point uj = of the integral fl4.36|) yields contributions that 
vanish for long t faster than any negative power of t. Just as in the transverse case, this 
result is a consequence of the fact that to this order in the HTL resummation, the spectral 
density is regular (no branch points) near uj = 0. 

To summarize, we gather the final results for the asymptotic real-time evolution of the 
transverse and longitudinal non-equilibrium expectation value in the linear approximation 

• Transverse (no external source): 



A Ti (k,t) = A Ti {k,0) 
+ A Ti (k,0) 



Z.[T]co sM -^^Uo(i 



(4.36) 



The sum rule ( [4.171 ) implies that the coherent field configuration relaxes to an asymp- 
totic amplitude which is smaller than the initial, and the ratio of the final to the initial 
amplitude is completely determined by the thermal wave function renormalization. 
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Longitudinal (impulsive external source): 

"asymp 

(k,t) 



Mk,t) = -Z l [T]^^ + aT t 



UJ 



l + O 



(4.37) 



The transverse and longitudinal wave function renormalizations Z*[T] , Z [T] are given 
by Eqs. ( ^T5|) and ( ^28|) respectively and a c ^ ymp {k,t) is given by Eqs. ( [PEIJOSD 
while the positions of the poles uj p , ojq can be obtained by the solving Eqs. ( |4.11| ) and 
( |4.24|) respectively. 

In Summary: we find that the long time dynamics is dominated by the Landau 
damping thresholds at u = ±k, not by the uo « region of the spectral density. The 
early time dynamics is determined by moments of the total spectral density. Long time 
power law tails had recently found in the current correlators in reference p8[ . 



V. LANGEVIN DESCRIPTION FOR THE SOFT MODES AND 
FLUCTUATION-DISSIPATION RELATION 

A. Langevin Description 

So far we have studied the damping of soft modes in the plasma by the hard particles 
from the microscopic point of view. In this section we provide a stochastic description of the 
relaxation of gauge fields via a semiclassical Langevin equation with a Markovian damping 
kernel and a gaussian white noise. 

A semiclassical description treats the hard modes as a 'bath' and the soft modes as the 
'system'. The bath degrees of freedom are integrated out, their main effect being encoded in 
a dissipative kernel and a stochastic noise inhomogeneity in the resulting Langevin equation. 
The dissipative kernel is related to the stochastic correlation function of the noise via a 
generalized fluctuation-dissipation relation. Physically the stochasticity arises because the 
hard scales which are integrated out in the HTL scheme and are responsible for Landau 
damping will also provide random kicks to the soft degrees of freedom. 

This section is devoted to a microscopic derivation of the Langevin equation for the in- 
homogeneous gauge field configuration, to leading order in the hard thermal loop program. 
This is achieved explicitly by integrating out the hard modes which provide a natural realiza- 
tion of the 'bath' variables while the soft modes are to be treated as the 'system' variables. 
The methodology of this approach is based on the Feynman- Vernon influence functional |50 | 



which has already been used to describe dissipation and decoherence in quantum systems 
from a microscopic theory p~i] Fi1f^ • 



The non-equilibrium partition function for the full field theory is 
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Z = J X?4? v ® {±) v ® {±)] ex P i J d *x {^ $(+)t d ^ (+) 



(5.1) 



e4 +) • - e 2 4 +) • 4 +) $C+)t $(+)) _ [+ 



and the effective action for the 'system' (the soft photons) follows by performing the path 
integral over the 'bath' (the hard scalars) treating the 'system' degrees of freedom as back- 
ground fields. This means that all subsequent expectation values will be evaluated in the 
reduced density matrix which defines the effective field theory for the system degrees of 
freedom. It is then convenient to introduce the 'center of mass' and 'relative' coordinates 



4 l = 



(±> 



R 



Ti 



(5.2) 



In terms of these redefined fields the effective action can be obtained to one-loop order via 
a systematic loop-expansion of the reduced partition function: 



S eff [A Ti , R Tl ]= J d 4 x {-R Ti d 2 A Tl - 2e 2 R Tl A Tl 



(5.3) 



le 



+ J <fx'R Ti {x) R TJ [(j^Wj^ix')) + OVMjSV)) 



;(-)« 



;(-)/ 



;(-)« 



Rtj (x 



d*x>R Ti (x) \ (jn( x )3Tj'( x ')) ~ (j£V)ir7V)>l Q ^ ~ x 'o) A Ti(x') 



;(+)< 



;(+)/ 



;(-)< 



The transverse current jTi{x) was introduced in Eq. fl2.2| ) and the curent- current correlators 
can be calculated easily using the defining formulas for the free scalar propagators in Eq. (|2l 



2. 11|) . In order to make explicit the soft momentum scales of interest we perform a spatial 



Fourier transform, in terms of which the reduced effective action, including the influence 
functional of the hard modes is given by 



&eff [An, Rn 



d 3 k 



dtl-Rn(k,t) (^ + k 2 + 2e 2 



($t$)J A Ti (-k,t)+ (5.4) 



(2tt) 3 

- /* dt'R Ti (k, t)Vij(k] t, t')A Tj {-k, t')+ 

+i J dtRn$,tWij(k',t,lf)RTj(-k, if 

where XV,- and J\T%j will be shown to be the dissipation and noise kernels respectively. The 
dissipation kernel which is given by 



T>ij(k;t, i') = 
d 3 p 



(5.5) 



4ie 



2e 



(2tt) 3 
d 3 p 



p Tl PTj [G<(t, t')G< +p (t, - G>(t, t')G> +p (t, t')Mt - t') 
PnPTj {(1 + n p + n p+k ) sm[(u p + u k+p )(t - t'))+ 



(2-K) 3 UJ p UJ k+p 

+ (n p - n p+k ) sin[(u; fc+ p - u p )(t - t')]} 6(i - t') 
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clearly gives a causal contribution to the effective action as seen in Eq.( |5.4| ). The fact that it 
is real follows from the properties of the non-equilibrium Green's functions [Eqs. (|2.10| ) and 
which imply that 

G<(t,t')G< +p (t,t>) - G>(t,t')G> +p (t,t') ~ 2zlm [G<(t,t')G< +p (t,t') 

Furthermore the dissipation kernel is in fact precisely the one-loop self-energy which appears 
in the effective action equation of motion Eq. (|3.6| ). The noise kernel on the other hand is 
acausal and gives an imaginary contribution to the effective action 



Uijik-t.t') -- 
d 3 p 



(5.6) 



(2tt) 3 
d 3 p 



p Tl PTj [G<(t, t')G< +p (t, + G>(t, t')G> +p (t, 0] 

Pn PTj {(1 + n p + n p+k + 2n p n p+k ) cos[(cj p + uj k+p )(t - t')]+ 



2 J (2ir) 3 u p uj k+p 
+ (n p + n p+k + 2n p n p+k ) cos[(u k+p - u p )(t - t')]} 



Therefore, 



Afij oc Re 



d 3 pp Ti PTj G<(t : t')G< +p (t,t') 



(5.7) 



which means that the noise-noise correlator and the dissipation kernel are the real and imag- 
inary parts respectively of the same analytic function of (t — t'). Thus they are automatically 
related by dispersion relations which reveal the fluctuation-dissipation theorem within this 
context. 

The imaginary, non-local, acausal part of the effective action gives a contribution to the 



path-integral that may be rewritten in terms of a stochastic field as |p2 
Z = J VR Ti (k) VR* Ti (k) exp 



d 3 k 

dt dt'—^ R* Tl (k, t) Mij{k\ t, f) R Tj (k, t') 



(5.8) 



ex J VR Ti (k) VR* Tl (k) VUk) VQ{k) P[£] exp 



d 3 k 



where the noise has a Gaussian probability distribution 

P[i] = exp[- J dt dt' d(k,t) K^(k;t,t') ^(-k,t')) 



(5.9) 



with zero mean and non-Markovian correlations: 

«&(M)»=0 ; «^(k,t)^(k',t') »= (2ir) 3 5 {3) (k + k') Afij&t^t'). (5.10) 

where <<•■■>> is the stochastic average with the probability distribution P[£\. 

In the Hard Thermal Loop approximation the leading terms in the noise kernel can be 
calculated explicitly, giving 
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Afij{k;t,t') 



e 2 T 3 



sm[k(t - t')] cos[k(t - t')] 
k 3 (t - t') 3 P(t-t') 2 



(5.11) 



Clearly these are long range correlations which cannot be replaced by local delta function 
in any regime of approximation. In fact the Fourier transform of the kernel yields: 



d{t - t') Nij{k;t - t') e 



-iuj(t-t') 



e 2 T 3 



T 



(5.12) 
(5.13) 



where £j(ia> + + ) is the imaginary part of the self-energy given by Eq. (|4.10[) resposible of 
the Landau damping processes. 

This clearly shows Landau damping to be the origin of the noise correlation and also 
provides an explicit realization of the fluctuation-dissipation theorem. We draw attention 
to the factor of T ' ju) which arises from the high T limit of the Bose factor (1 + 2n UJ ) (see 
below) and leads to a noise-noise correlation oc e 2 T 3 . This result is in accord with the usual 
fluctuation-dissipation relation, in which the noise-noise correlation function has one more 
power of T than the dissipative contribution to the Langevin equation. 



The Langevin equation itself is obtained via the saddle point condition [32 



SS, 



5Ri{k,t) 



\R=0 



(5.14) 



leading to 



dt 2 



-, rt -. 

A Ti (k,t) + k 2 A Ti (k,t) + 2e 2 (&$)A Ti (k,t) - / dtfT>ij(k; t, t')A T j{k, t') =&(M') 

(5.15) 



B. Fluctuation-Dissipation relation: 

The general form of the fluctuation-dissipation relation can be established at this level by 
retaining the complete expressions for the dissipative and noise kernels without performing 
the HTL approximation. A standard analysis of the spectral representation P^| , |55| of the 
equilibrium correlators leads to the following result for the Fourier transform of the noise- 
noise correlation function 



Afijik, 



- Im 

4 



coth 



~2~ 



(5.16) 



Finally taking temporal Fourier transforms of both sides of the Langevin equation ( |5.15 ) 
and averaging over the noise with the Gaussian distribution function we find 
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<< ATi(k,co)Ari(—k, —J) »= \A$i(k, tu)\ 2 + p(k, u) coth 



(3u 



I Im 






4 (u 2 -cu 2 k -Re 


'^f(k,u\ 


) 


2 + (im 




r 



(5.17) 



where the tildes stand for the Fourier transform and A^^k, t) is a solution of the homoge- 
neous equation, which is precisely given by Eq. ( |4.36| ) (assuming A^(k, t = 0) = which can 
be relaxed with the proper generalization). The double brackets stand for averages over the 
noise with the probability distribution (5T§)- fl5TT0D . This is the form of the usual fluctuation- 
dissipation relation, which we obtained consistently by integrating out the hard modes and 
deriving the influence functional [H^fyJ for the transverse components of the gauge invariant 
fields. The semiclassical Langevin equation is useful in order to obtain semiclassical corre- 
lation functions by averaging the solution of a partial differential equation over a stochastic 
gaussian noise. The question arises: to which correlation function of the microscopic theory 
are these stochastic averages related? The answer to this question is found by writing a 
spectral representation of the equilibrium correlator {A i ^(t)A i _?(£')) in which the brackets 
stand for averages in the equilibrium density matrix. A straightforward but tedious exercise 
following the steps described in [P5| , p^1 reveals that this relation is given by 

« A Ti {k,uj)ATi{-k, -u) »= - {(Ai(k,u)M-k, -<*>)> + (A{-k, -oj^ftw))} (5.18) 

which again is a result known to be a consequence of the fluctuation-dissipation relation in 
simple systems. The correlation function (|5.18| ) has a finite, non-trivial classical limit and 
agrees with the one proposed to be studied within the context of classical field theory in [E3| . 



The real-time analysis presented here agrees with the general picture discussed in refs. 
29| . Moreover, our analysis reveals the precise range of the kernels. 
Clearly the situation will be more complicated in QCD where a separation between hard 
and soft degrees of freedom must be implemented in order to obtain the influence functional 
for the soft degrees of freedom. However, the procedure detailed in this section can be carried 
out consistently once this separation is introduced. 

Of course the main rationale for obtaining a Langevin equation is to provide a semi- 
classical scheme to implement the calculation of correlation functions from the solutions of 
stochastic differential equations. However, we note that unless a successful scheme to deal 
with the non-Markovian kernels is implemented the advantages of a Langevin description 
are at best formal. Ignoring the non-localities of the dissipative and noise kernels will clearly 
miss the important physics associated with Landau damping. A naive Markovian approx- 
imation is not only uncontrolled and unwarranted but clearly very untrustworthy in view 
of the fact that the relevant kernels are all long-ranged and it is precisely this long-ranged 
nature of these kernels which is responsible for the important dissipative effects of Landau 
damping. The importance and difficulties of keeping these non-localities in a classical lattice 
description has been recognized in [j33| -|35|l- 
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VI. KINETICS OF LANDAU DAMPING 



The real-time formulation of non-equilibrium quantum field theory allows us to obtain the 
corresponding kinetic equations for the relaxation of the occupation number or population 
of quanta. In particular our goal is to obtain the kinetic equation for the relaxation of the 
expectation value of the number of soft quanta. In keeping with the focus of this article we 
will only consider the population relaxation to lowest order in the HTL approximation and 
concentrate mainly on the understanding of relaxation via Landau damping. 

Kinetic approaches towards describing transport phenomena and relaxational dynamics 
typically require a wide separation between microscopic time and length scales, namely the 
thermal (or Compton) wavelength (mean separation of particles) and the relaxation scales 
(mean free path and relaxation time). This approach which ultimately leads to the Boltz- 
mann transport equations involves the identification of slow and fast variables which justifies 
a gradient expansion. This is a coarse graining procedure that averages over microscopic time 
scales and leads to irreversible time evolution. 

In the collisional approach to Boltzmann kinetics only the distribution-changing processes 
that conserve energy and momentum are considered, and these are weighted by the corre- 
sponding Bose/Fermi statistical factors. Off-shell processes that occur on time scales shorter 
than the relaxation scale are not included; therefore only processes with asymptotically on- 
shell final states are accounted for in this description. Landau damping processes which 
contribute via thermal loops to forward scattering are not included in the typical Boltzmann 
equation. Thus we anticipate that the naive Boltzmann approach will yield no population 
relaxation via Landau damping. However this is conceptually inconsistent because we have 
learned in the previous sections that an initial coherent configuration will relax to an am- 
plitude smaller than the initial by these processes and we would expect such a relaxation 
to contribute to a depletion of the number of quanta of the initial state. The resolution of 
this inconsistency requires one to go beyond a simple Boltzmann approach and to include 
off-shell processes in the kinetic description. This is the focus of this section. 

The contribution of off-shell processes to the non- equilibrium evolution of particle distri- 
butions is ignored in most approaches towards kinetics. Recently an approach to kinetics 
that incorporates off-shell effects has been proposed to describe processes in which relax- 
ation competes with other fast scales pOHfl| , in particular near phase transitions [ 40[ . The 



importance of off-shell contributions to the evolution of particle distributions has also been 
recognized within the context of fast kinetics in semiconductors ||42|| . 



This section is devoted to a study of kinetics as an initial value problem pR3j-f42"l in order 
to reveal the role played by off-shell processes in transport phenomena and relaxation in the 
medium. 

We derive a kinetic equation that takes into account microscopic time scales in the 
theory from first principles allowing us to analyze clearly the effect of off-shell non-collisional 
Landau damping processes on the evolution of the photon distribution. This framework will 
allow us to make contact with the relaxation of a coherent initial configuration studied in 
the previous sections. We want to study both the relaxation of an initial distribution of 
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asymptotic photons with free dispersion relation Q(k) = k, as well as for the quasiparticles 
with dispersion relation Q(k) = u p (k) with u p (k) being the solution of the dispersion equation 
( [4.1 1| ) i.e. the 'true', in-medium pole. We will distinguish between these two physically 
different cases and address them separately. 

This approach begins by defining a suitable number operator 



|D|]. In the case of asymp- 



totic photons this is the usual number operator in terms of the canonical field and momenta 
which is given by the energy per momentum k divided by the frequency. In the case of 
quasiparticles the energy stored in the plasma has two components: the free field part plus 
the response from the medium. From classical electromagnetism of polarizable media ]54| 
in linear response, the two contributions lead to an energy density in the medium that is 
quadratic in terms of the electric and magnetic fields, each term however, multiplied by a 
coefficient that involves derivatives of the dielectric and permeability tensors with respect to 
frequency ||54|| . The contribution from the plasma collective modes is obtained by evaluating 



these coefficients at the plasma frequency. Using Kramers Kronig dispersion relations [153 



these coefficients are related to the residue of the dielectric constant at the plasma poles. This 
relation has been formalized at the field theoretical level by Migdal in his pioneering work 
on collective modes in medium where he developed the quantization procedure in medium 
in terms of quasiparticle operators ||55| . This field theoretical treatment automatically leads 
to the identification of these coefficients with the residues of the retarded propagators at 
the plasmon poles, i.e. the wave function renormalization. Migdal obtains in this manner 
the energy density corresponding to on-shell collective modes in terms of the quasiparticle 
operators |55[ . 



More recently the energy density of the plasma including the polarization effects has been 
obtained in terms of operators that create and destroy collective modes in the plasma by 
Blaizot and Iancu [p6 |. The results of these authors is consistent with the collective mode 
quantization and the energy density obtained by Migdal |SJ| and with the results of classical 
polarizable media [54]. Blaizot and Iancu |56| use the collective mode decomposition of the 
transverse gauge field 



\ 2u p (k) 



£ [e x (k)a x (k)e~^ k)t + e x {-k)a{{-k)e l ^ k)t 



(6.1) 



A=l,2 



Where the operators <x\{k) ; a x (k) have a retarded propagator with unit residue at the 
plasmon pole and Z\ is the wave function renormalization. 

Migdal |5_5] and Blaizot and Iancu [5B] prove that the energy density associated with the 
collective modes in the medium can be written as 



£{k)=u p (k) a\(k)a x (k) 



(6.2) 



A=l,2 



This result is the same as that obtained in the classical theory of polarizable media when 



the electric and magnetic fields (transverse) are written in terms of collective modes [54 



Thus following Migdal |55[ and Blaizot and Iancu |56j we introduce the Heisenberg num- 
ber operator of on-shell collective modes 
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A T {k) ■ A T (-k) 



Q 2 A T (k) 



M-k) 



2Z{ 



(6.3) 



where for asymptotic photons Q k = k 
contribution) and for collective modes Qk 



1 (we neglect here the zero temperature 
uj p {k) ; Z\ = Z t \T\ withcap(fc) being the plasmon 
pole and Z l \T\ the wave function renormalization in the HTL limit. This formulation now 
permits to treat the collective modes much in the same manner as the usual renormalized in 
and out fields in S-matrix theory, i.e. by rewriting the action in terms of the renormalized 
fields and introducing counterterms that reflect the true position of the pole and residue. 

It turns out to be easier to work with the Heisenberg operator N k rather than the number 
operator and is obtained using the Heisenberg operator equations which are easily seen to 
be: 

JTt(M) =fir(-M) 

U T (-k, t) = -{k 2 + 2e 2 ($ t $) + 5n 2 k )A T (k, t) - j T {k, t) . 

where the counterterm accounts for the definition of the number of quasiparticles and Hr(/c) 
represents the canonical momentum conjugate to the transverse electromagnetic field Ax{k). 

We now consider an initial state described by a density matrix for which the expectation 
value of the above number operator for on-shell collective modes is non- vanishing. 

Using these equations the expectation value of the Heisenberg rate operator in the initial 
density matrix is obtained in the following form which is rather convenient for subsequent 
calculations: 



(N k )(t) = N k (t) 



1 







t=t ,: 



2e 2 ($ t $) + 5Q 2 k d_ 
' W 1 



(j+(k,t) ■ A^{-k,t")) 
(A+(k,t) ■ A T {-k,t") + (A+{k,t") ■ Ar(-k,t)) 



t=t" 



(6.4) 



This formulation has been previously applied to the study of photon production in a strongly 
out of equilibrium phase transition [|40[j . 

The expectation values are calculated by inserting the operators into the closed time path 
integral and expanding in powers of a. Since 5Vt\ is of order a the second term in ( |6.4j ) is 
calculated as a tadpole in free field theory and it vanishes identically. Let us consider the 
case in which the initial density matrix at time t is diagonal in the basis of eigenstates of 
the number operator, and evolves subsequently with the interaction Hamiltonian. To lowest 



order in a we find the expectation value of the rate to be pi, 



N k (t) 



alt' / dp [p 2 - (p ■ kf 



(6.5) 



x 



G>(t,t')G> +p (t,t')g<(t',t) -G<(t,t')G< +p (t,t')g>(t',t)] Q(t 



t') 



28 



The theta function ensures that this expression is causal. The Green's functions for the 
scalars and the photons can be read off from Eqs. (|2.10|) , (|2.11| ),( [2.12j ) and (|2.13|) but with the 



frequency Q k replacing the bare frequency. Finally the rate can be written as pEOfl : 

N dt)= TTT^I [ ^P 2 si^6 fdr{ 

COs[(oJ p + UJk+p + Q k )(t - r)] [(1 + N k {t )){l + n p ){l + n k+p ) - N k (t )n p n k+P ] + 
cos[(uj p + u k+p - Q k )(t - r)] [(1 + N k (t ))n p n k+P - N k (t )(l + n p )(l + n k+p )} + 
cos[(^ p - u k+p + Q k )(t - r)] [(1 + N k (t ))(l + n p )n k+p - N k {t )n p (l + n k+p )} + 
cos[(^ p - iu k+p - Q k )(t - t)] [(1 + N k (t Q ))(l + n k+p )n p - N k (t )n k+P (l + rip)]} (6.6) 

We note that the expression above depends on the occupation number of the gauge field 
at the initial time to only - this is obviously a consequence of the fact that perturbation 
theory at lowest order, neglects the change in occupation number. Recently we have pro- 
posed |^0| a Dyson-like resummation of the perturbative expansion that includes off-shell 



effects in the relaxation of the population. This resummation scheme is obtained by the re- 
placement iVfc(to) — ► N k (r) in (|6.6|) resulting in a non-Markovian description. The resulting 
kinetic equation with memory is akin to that obtained via the generalized Kadanoff-Baym 



approximation [42] in non-relativistic many body systems. This approximation has been 
recently shown to coincide with the exact result in the weak coupling limit in a solvable 
model of relaxation [0,0 and will be shown below to imply a Dyson-resummation of the 



perturbative series. The kinetic Eq. (|6.6|) has an obvious interpretation in terms of gain 



minus loss processes ||40|| , but the retarded time integrals and the cosine functions replace 
the more familiar energy conserving delta functions. Taking the occupation number outside 
the integral and integrating to large times, thereby replacing the cosines by delta functions 
as in a Boltzmann description would lead to a vanishing right hand side since none of the 
resulting energy conserving delta functions can be satisfied. However, the non-Markovian 



kinetic equation ( |6.6j) will lead to non-trivial relaxational dynamics that will be studied in 
detail below. 

Let us consider the situation in which the initial state has been prepared far in the past, 
i.e. to - ► — oo. 

An equilibrium solution is simply 

N^ q = constant (6.7) 

where the constant is arbitrary because none of the resulting energy- conserving delta func- 
tions can be satisfied for the values of Q k either corresponding to the bare frequencies or the 
quasiparticle poles. This is a consequence of the off-shell processes, a detailed understanding 
of this feature will be provided elsewhere [58.^9 . 



Let us now consider departures from this equilibrium solution and study the relaxation 
of a disturbance in the distribution function introduced in the system at t = so that 

N k (t = 0) = N e k q + 5N k (0). (6.8) 
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Denoting the particle distribution for t > by 

N k (t > 0) = N? + 6N k (t) (6.9) 

we obtain a rate equation for 5N k (t) which is now the same as in the previous step except 
that the time integrals stretch from to t instead of to —* — oo to t 

d 6^ f d^v ft 

T + 6N k(t) = 1fi 30 ~ t / —P 2 sin2 6 dr{(l+n p + n p+k ) 

dt lo7r i SZ fe .Z| J u p u k+ p Jo 

[cos[(u p + Lu k+P + Q k )(t - r)] - cos[(u p + u k+p - Q k )(t - r)]] 

(n p+k - n p ) [cos[(cj p - uj k+p + n k )(t - r)] - cos[(cj p - uj k+p - Q k )(t - r)]]} 5N k (r) . (6.10) 

Terms independent of N k vanish identically since the time integrals for those terms yield 
delta functions which are never satisfied. 

Although this non-Markovian but linear equation will be solved exactly by Laplace trans- 
form below, it is illuminating to compare the different approximations that are obtained un- 
der the assumption that the relaxation time scale for SN k is much longer than the time scale 
of the non-local kernels. Under this assumption, which will be analyzed below, 5N k (r) can 
be replaced by SN k (t) and taken outside of the integral leading to a Markovian description. 
A further approximation, taking the upper limit of the remaining integral to t — > oo leads 
to the familiar Boltzmann equation, thus the two approximations to be compared with the 
'exact' solution are the following 

• Markovian: 

pN"(t) = -T(t) 5N k M (t) (6.11) 

r ^ = -J7r4v^t I ~^^P 2 sin 2 9 f dr {(1 + n p + n p+k ) 
[cos[(^ p + u) k+p + Q k )(t - r)] - cos[(^ p + uj k+p - Q k )(t - r)]] 

(n p+k - n p ) [cos[(ujp - u k+p + Q k )(t - r)] - cos[(cj p - u k+p - Q k )(t - r)]]} , (6.12) 
with solution 

5N™(t) = e -Jo tr W 6Nj?(0) . (6.13) 



• Boltzmann: 



^8N*(t) = -r(oo) 6N*(t) (6.14) 
8N?(t) = e- r(oo) * 5N k B (0) . (6.15) 
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Taking the limit t — > oo in T(t), the cosines become energy conserving delta functions 
and comparing with the expression for the transverse self energy given by Eq.( f4.7|) , it is 
straightforward to see that in the Boltzmann approximation we obtain 



dt 



6N*(t) 1 



f2 fc 



{t) = 



(6.16) 



with Sj(fifc) being the imaginary part of the transverse self-energy evaluated at Qk- This 
result is the familiar relationship between the relaxation rate of the particle distribution 
r(oo) and the damping rate, which is determined by the imaginary part of the self-energy 
on-shell and vanishes in the present case because the damping processes are off-shell. 

The solution of the non-Markovian equation ( |6. 10 ) is obtained by Laplace transform and 
given in general by 



rc+ioo A a 

5N k (t) = 5N k (0) / — 



S k {s) 



(6.17) 



where c is a real constant chosen so that the contour is to the right of the singularities of 
the integrand and Skis) is the Laplace transform of the non-local kernel given by 



Skis] 



d3 P 2 ■ 2 a 
p sin v 



+ 



167T 3 fi fc J UJpUJk+p 

s 

S 2 + iuJ p + UJk+p - fife) 

s 



s 2 + [uj p + u k+p + fifc) 2 
(1 + n p + n k+p ) 



_S 2 + - UJk+p + fifc) 2 S 2 + (U p - UJk+p - flkf 



i n p+k — n p) 



(6.18) 



The expression ( |6.17|) clearly shows that the non-Markovian rate equation ( |6.10|) implies 
a Dyson-like resummation of the perturbative series as anticipated before. 

Comparing the non-Markovian rate equation ( |6.10|) to the Markovian approximation 
( |6.11f) one can clearly see that the Markovian approximation averages over the time scales 
of the kernel, whereas the non-Markovian equation includes the contribution of coherent 
processes throughout the history of the kernel. If the range of the kernel was indeed shorter 
than the relaxation time of the population, the real-time solutions of both equations will 
differ by terms of the order of the ratio of the time scale of the kernel to the relaxation 
time scale of the distribution. However, in situations in which the kernel is long-ranged as is 
the case under consideration, the non-Markovian expression allows the inclusion of coherent 
effects in the relaxation. 

The dominant contribution to S{s) in the HTL limit arises from the Landau damping 
term leading to the simplified expression 



Skis) 



e 2 T 2 k s 



x(l 



X 



1 s z 



[kx — flkY 



dx 



(6.19) 
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More explicitly, 

e 2 T 2 f A s(k 2 + s 2 - 3ft 2 ) , s 2 + (k - ft,) 2 

+ — 2^ log ( a -^ + ng ) • (6 - 20) 

It is easy to show that this form of Sk{s) is in fact related to the self-energy in the HTL 
limit by the following illuminating equation: 

„ . , E*(s -in h ) - E*(s + iQ k ) 

Sk(s) = ^ (6.21) 

The real-time solutions to Eqs. (|6~l0| ), (|6U7|) and (|6TTg ) in the HTL limit with ( gig ) 
will now be given for the cases of quasiparticle and bare particle respectively. This analysis 
will reveal the range of validity of the assumptions leading to the Markovian approximation. 



A. Quasiparticle: Qj- = Wp(k) 



When Qk is chosen to be the quasi-particle pole, uj p (k) > k we note that the integrand of 
( |6.17|) has a single isolated pole at s — 0. Indeed, the limit of s — > in Sk(s) would lead to 
delta functions in Eq. fl6.19|) which, however, cannot be satisfied for Qk = &p(k). Therefore, 
<Sfe(0) = and s = is an isolated single pole and completely determines the asymptotic 
limit of the real-time solution. This can also be seen by looking at Eq. ( p.21] ), from which 
the analytic structure is explicit. Clearly, Sk(s) vanishes at s = because £*(— iuj p {k)) = 
Y?(iu p {k)). Furthermore from the known singularity structure of E*(s) one concludes that 
<Sfe(s) must have branch cuts for i{ui p {k) — k) < s < i(u p (k) + k) and —i{uj p {k) — k) < s < 
-i(w p (k) + k). Using Eq. flPID , it is now a straightforward exercise to see that the residue 
at this pole at s = is given by (1 -2dY> t (iVL)/dVL 2 y 1 which to this order is w Z*[T] 2 where 
Z l [T] the (transverse) wave function renormalization given in Eq. ( [4.151) . 

Again the long time behavior is completely dominated by the end points of the cut, 
leading to the asymptotic result 



5N k (t) t = c "5N k (0){Z t {T] + 



e 2 T 2 7r 2 



YlZ l \r\ku) p {k)t 2 
cos (u p (k) — k)t 



(ou p (k) - k) 2 (l + 



e 2 T 2 _D_ 
3k 2 



cos (u} p (k) + k)t 



,(*)+*)» (1 + ^) 



[UJ 



(6.22) 



where 
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We clearly see that asymptotically the population has relaxed to a smaller value and the 
ratio of the asymptotic to the initial population is determined by the square of the thermal 
wave-function renormalization. This is in agreement with the analysis of the relaxation of the 
expectation value of the field - since N k oc A\ it is expected that the ratio of the asymptotic 
value of the quasi-particle population to the initial value be proportional to the square of 
the same relation for the expectation value of the field. 

Therefore the relaxation of the quasiparticle number has its origin in Landau damping, 



this is consistent with the results of Blaizot and Iancu [ 56| who proved that the time derivative 
of the total energy is related to Landau damping. Since the quasiparticle number is related 
to the energy of the collective modes the relaxation of the quasiparticle number is directly 
related to Landau damping. 

This intuition is also borne out in the usual Boltzmann approach wherein the relaxation 
rate of the distribution function (in the relaxation time approximation) is twice the damping 
rate for the quasiparticle. 

This relaxation has a simple interpretation. Consider the case of a physical electron with 
interpolating operators defined to create single electrons with the physical mass and unit 
amplitude out of the in or out vacuum states. The asymptotic correlation function in real 
time of these interpolating operators has the oscillatory parts corresponding to the physical 
pole (with unit residue), but there are power law corrections arising from the overlap of 
the states created by these operators with the multiparticle continuum [pTfl . Whereas at 
zero temperature the multiparticle continuum is beyond the two particle threshold, at finite 
temperature and in the case under consideration, the leading contribution is obtained from 
Landau damping corresponding to intermediate states with space-like momenta. 

For the Markovian approximation (|6.13| ) we find for t — > oo 
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(6.23) 



For e 2 T 2 /k 2 << 1 we see that to lowest order in e 2 T 2 /k 2 the perturbative expansion of the 
Markovian solution coincides with the solution of the non-Markovian equation. However for 
soft momenta such an expansion is not valid and the validity of the Markovian approximation 
must be questioned. 



B. Failure of the Markovian and Boltzmann approximation: 

To assess whether the Markovian and Boltzmann approximations will be reliable we must 
understand the different time scales, in particular the range of the kernel. 

In the Hard Thermal approximation we find that the non-Markovian kinetic equation 
( p.l(J| ) reduces to, 



^-5N k (t) ~ t dx (1 - x 2 )x f dr sm{kx(t - r)) sm[Q k {t - r)} SN h {r). (6.24) 

dt l^i'fc J -i Jo 
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The integral over the variable x inside the kernel can be performed and we find that the kernel 
falls off as l/(t — t) 2 + • • •. Then if eT/k << 1 the relaxation time scale of the population is 
longer than the range of the kernel and the Markovian approximation is warranted. In this 
case the discrepancies between the non-Markovian and Markovian results are perturbatively 
small. On the other hand, for soft scales the relaxation time scales become comparable to 
the time scale of the kernel and a Markovian approximation is certainly unjustified. The 
non-Markovian equation for relaxation includes the coherent effects on similar time scales 
that are averaged out (coarse-grained) in the Markovian approximation. 



C. Bare particle, or hard quasiparticle: 



For the case of the bare particle, the dispersion relation is simply £l k = k - this is also 
the case for the large k limit of the quasiparticle dispersion relation pl| , |32| , |37 |. In this case, 
the quantity S k (s)/s in ( p,18| ) has a logarithmic singularity as s — > 0, because the position 
of the putative pole fl k has moved to the tip of the cut and there is a pinching singularity. 
There is no longer a pole at s = in the Laplace transform [s — S k (s)}~ 1 , rather it diverges 
as (sins) -1 as s — - > 0. This logarithmic divergence arising from the pinching singularity is 
very similar to that recently studied within the context of hard fermions |HIJ . 

In the Markovian approximation ( |6.11| ) and in the Hard Thermal Loop limit, the rate 
equation becomes 
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(6.25) 



which for long times yields a power law with an anomalous exponent: 

2 2 

6N k (t) ~ SN k (Q)(kt)-^. 



(6.26) 



An anomalous exponent somewhat similar to this one has been found in [60] in the case of a 



hard fermion and it has the same origin, i.e. a pinching infrared singularity. The expression 
( |6.25| ) reveals that the kernel is long ranged, falling off with an inverse power of time in this 
case. Therefore a Markovian approximation will be justified when e 2 T 2 /k 2 « 1 because 
only in this weak coupling limit is the population relaxation slower than the fall off of the 
kernel. 

The solution of the non-Markovian equation ( |6.10| , |6.18| ) in the HTL limit (only the 
Landau damping contribution is considered) is obtained again by inversion of the Laplace 
transform. We find that the long time behaviour is given by the end-points. In this case the 
uj = end-point dominates yielding 



8N k (t)=5N k (0) 
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(6.27) 



where e = 2.718 ... is the base of the natural logarithms. 

It is illuminating to point out that as ~~ > & the logarithmic singularities in the real 
part of the transverse self energy imply that (for fixed finite k) 

^^Uk->oo =► Z*[T]^0 (6.28) 

As the pole approaches the tip of the cut (for finite k) the residue becomes smaller until it 
vanishes exactly when the pole merges with the continuum. It is remarkable that in this limit 
distortions of the distribution function relax completely and vanish asymptotically. Although 
for the case of a hard quasiparticle with k » eT the dispersion relation approaches that 
of a free particle, Z f [T] — > 1 and physics is perturbative, there is, however, slow relaxation. 
We note that none of these effects can be captured by a simple Boltzmann approach since 
all of these phenomena are associated with off-shell effects. 

In the case of the distribution function for bare particles, we can interpret this anomalous 
relaxation as the dressing effect from the medium, i.e. at long times the bare particles are 
completely dressed by the medium and disappear from the spectrum. 

The interpretation is different for the case of a hard quasiparticle, in which case the free 
dispersion relation is obtained in the limit k » eT. In this limit the effective coupling 
eT/k << 1 and the relaxation is slow. This result could be important in understanding 
the relaxation of a distribution of photons produced in Bremmstrahlung processes at high 
energy in the quark gluon plasma. A full study of the relaxation of hard quasiparticles is 
beyond the realm of this article and will be studied in detail elsewhere ||59||. 



VII. DISCUSSION AND CONCLUSIONS 

Our goal in this article was to provide a detailed analysis of the real-time relaxation of soft 
gauge invariant non-equilibrium expectation value through the off-shell process of Landau 
damping. These determine the leading contributions to the thermal propagators in the Hard 
Thermal Loop limit and are the dominant contributions to the long-time asymptotics. The 
off-shell nature of these processes determine the non-Markovian nature of the relaxation 
phenomena associated with them. We focussed our study on the leading HTL contributions 
to the relaxation of gauge invariant transverse and longitudinal non-equilibrium expectation 
value in scalar electrodynamics. These results will also apply to fermion electrodynamics 
and non-abelian theories since the structure of the retarded Green's function is the same 
up to this order in the HTL expansion. After providing an elegant re-derivation of the well 
known HTL effective action using the tools of non-equilibrium quantum field theory in the 
linear amplitude approximation f36|, we moved on to the main goals of this article: 
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• i: to study in detail the relaxation in real time, in the linear approximation (small am- 
plitude) of the transverse and longitudinal gauge invariant non-equilibrium expectation 
value. The off-shell process of Landau damping results in power law relaxation of the 
transverse and logarithmic relaxation of the longitudinal (plasmon or charge-density) 
excitations. Both types of non-equilibrium expectation value relax to an asymptotic 
amplitude that depends on the thermal wave-function renormalization, which is com- 
pletely determined by Landau damping in the HTL limit. One of the main conclusions 
of this detailed analysis is that the relaxational dynamics asymptotically at long times 
and to leading order in HTL resummation is completely determined by the behavior 
of the spectral density near the Landau damping thresholds at uj = ±k (a branch cut 
singularity), the contribution from the region uj « are regular (no branch point sin- 
gularities) and therefore lead to subleading corrections to dynamics at long times and 
high temperature. This result is for both longitudinal and transverse non-equilibrium 
expectation value and is confirmed by an exhaustive analytic and numerical study. 
The short time evolution of the non-equilibrium expectation value is determined by 
moments of the total spectral density. Therefore a complete understanding of the the 
global analytic structure of the retarded propagator, in particular the complete cut 
contribution from Landau damping processes is required. 

This is special to the HTL resummation at one-loop: at higher orders, a branch point 
could develop at uj = 0. Such a branch point would produce pure power like tails, with 
no oscillations, and so dominate at large times. 

We restricted ourselves in this paper to small amplitude non-equilibrium expectation 
value so that we were confined to the linear regime. New phenomena beyond the HTL 
scheme are expected in the non-linear amplitude regime. Such regimes can be studied 
within self-consistent Hartree-type approximations in the out of equilibrium framework 

0- 

• ii: 

We have obtained the influence functional (non-equilibrium effective action) for the 
soft gauge invariant degrees of freedom by integrating out the hard degrees of freedom 
to leading order in the HTL approximation. This allowed us to obtain the Langevin 
equation for the soft degrees of freedom to leading order in HTL and to provide a 
microscopic ab initio calculation of the dissipative and noise kernels in the HTL limit. 
Both kernels display the non-localities associated with Landau damping and we find 
that there is no region of time scales in which a Markovian approximation describes 
the dynamics correctly. As a byproduct we obtained the fluctuation-dissipation rela- 
tion and recognized the correlation function that emerges in the classical limit. We 
established in detail that a Markovian description of relaxation of transverse or lon- 
gitudinal non-equilibrium expectation value is unwarranted to this order in the HTL 
resummation. 
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iii: Having understood the relaxation of coherent non-equilibrium expectation value 
through off-shell effects of Landau damping we asked how these processes can be in- 
coporated in the relaxation of the distribution function for the transverse fields. A 
Boltzmann approach would yield no relaxation to lowest order in the HTL because 
there is no imaginary part on-shell for the transverse or longitudinal quasiparticles. 
Therefore a kinetic description of the relaxation of the distribution function must nec- 
essarily go beyond a Boltzmann collision approximation. 

We provided a novel description of the kinetics of relaxation of the distribution function 
for transverse degrees of freedom that includes off-shell effects and goes far beyond the 
usual Boltzmann approach. We have compared the relaxation obtained from this non- 
Markovian kinetic equation to that described by the usual Boltzmann (which yields 
a trivial result) and a Markovian version that includes coarse-grained details of the 
off-shell processes. We have found that the distribution function for soft quasiparticles 
relaxes as a power (1/t 2 ), and found an unusual dressing dynamics for bare particles. 
This kinetic approach also reveals unusual logarithmic real time relaxation for hard 
quasiparticles resulting from infrared pinching singularities similar to those found in 
the case of hard fermions [6^ . 



The body of these results reveals new and unusual features of relaxation of soft degrees 
of freedom in gauge theories. These will obviously have to be taken seriously into account 
in a full description of relaxational processes in the QGP and should also be important to 
clarify better the role played by damping in the sphaleron rate in the symmetric phase. 

In higher order corrections there will arise contributions from collisional processes that 
provide an imaginary part on-shell. At next order for example both Compton scattering and 
pair-annihilation will contribute to collisional relaxation and will provide a collisional width 
both to the transverse and longitudinal (plasmon) degrees of freedom. The imaginary part 
of the self-energy on-shell is typically associated with a damping rate and associated with 
an exponential decay of the amplitude. Excepting intermediate time scales, this exponential 
relaxation, however is not a proper description either at early or long times where power 
laws dominate the dynamics ||36|| . The relaxation associated with Landau damping at lowest 



order, will have to be balanced with the next order corrections which yield an approximate 
exponential relaxation and the resulting dynamics will depend on the details of the com- 
petition between these different processes, both Landau damping and collisional. Which 
process dominates will depend on the particular time scale of interest and the time scales for 
competition between the two different type of phenomena will depend on the details of the 
perturbative contributions. The kinetic approach introduced here could also prove useful to 
study the energy loss of quarks and leptons via off-shell processes in the QGP. 

We plan to address this competition between Landau damping and collisional phenomena, 
along with an extension of the treatment presented in this article to leptons and a more 
detailed study of non-Markovian kinetics in future work. 
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FIG. 5. t x (ln[t]) 2 - 5 x A^ ut (k,t) vs. t (in units of 1/k) for m 2 D /k 2 = 2. 
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